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" . Abstract 

m . 

, We discuss the stability of (charged) static black holes in higher- 

' dimensional spacetimes with and without cosmological constant by us- 

ing gauge- invariant master equations of the Schrodinger equation type 
for black hole perturbations derived by the authors recently. In particu- 
lar, we show that the stability of higher-dimensional Schwarzschild black 
holes can be proved with the help of a technique called S-deformation of 
the master equations. We also point out that higher-dimensional static 
black holes might be unstable only against scalar-type perturbations in 
the neutral case and in the charged case with spherically symmetric or 



' flat horizons. 

1 Introduction 



The perturbation analysis of 4-dimensional black holes has a long history 
and has provided valuable information on various problems such as the 
fate of gravitational collapse, the stability and uniqueness of black holes, 
cosmic censorship, and gravitational wave emissions. From a technical 
point of view, the key point in this perturbation analysis was the fact that 
the perturbation equations can be reduced to a single 2nd-order ordinary 
differential equation (ODE) for a master variable $ of the form 

-0+[nO--1* = 5. (1) 



*A talk given at the Seventh Hungarian Relativity Workshop, 10—15 Aug. 2003, 
Sarospatak, Hungary. 
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Recently, starting from a gauge-invariant formalism for perturbations 
of higher-dimensional spacetimes with spatial symmetry we showed 
that similar reduction to single master equations is also possible for pertur- 
bations of generalised static black holes with or without charge in higher 
dimensions [3 01 • Further, using these master equations, we analysed the 
stability of such black holes |1] 01 • 

In this talk, we explain the basic idea of this formulation and describe 
the results of the stability analysis of (charged) static black holes in higher 
dimensions. 

2 Tensorial Decomposition of Perturba- 
tions 

2.1 Unperturbed Background 

We assume that the unperturbed background spacetime is locally written 
M ~ JV X IC^ and its metric has the form 

ds^ = 9ab{y)dy°-dy'' + r^{y)dal, (2) 

where dcrj; — ')ij{z)dz^dz^ is a metric of an n-dimensional complete Ein- 
stein space K,^ with Rij = (n — l)Kjij {K — 0, ±1), and gab is a static 
metric of the 2-dimensional space 7V^ expressed as 

dr'^ 

Qabdy'^dy^ = -f{r)dt^ + ——. (3) 
fir) 

In this and the next section, we assume that this metric is a solution to the 
vacuum Einstein equations with cosmological constant A — n{n + l)A/2. 
Hence, /(r) is expressed as|S] 

/M = ^-|^^-Ar^ (4) 

When the spacetime contains a regular black hole, the space IC" rep- 
resents a spatial section of the horizon. For K — 1, K," is always com- 
pact from Meyers' theorem and for n < 3, IC" is locally isomorphic 
to 5"*. The spacetime contains a regular black hole if AM^^'""^' < 
(n - l)/(n + In contrast, ior K = or -1, K." may not be 

compact, and the spacetime contains a regular black hole only for A < 0. 

2.2 Perturbations 

In terms of the perturbation variable -0^^ = hf_iu~hg^i, /2 with h^i, = (Jg^^, 
the perturbed vacuum Einstein equations are written as 

-VVm. - 2/?,,,,^^"" + 2V(„V"t/.,)^ - V" VVc/sffM. = 0. (5) 

In order to analyse the behaviour of perturbations using this equation, 
we must solve two problems. First, this equation is invariant under the 
gauge transformation 

x'^^x^'+e ^ = -VmC. - V.Cm- (6) 
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In order to extract the dynamics of the physical degrees of freedom, this 
gauge freedom should be eliminated. Second, @ is actually a set of 
coupled equations with (n+2)(n+3) /2 entries and is quite hard to analyse 
in general. 

For the special background ((2J, the second problem can be made 
tractable with the help of the following tensorial decomposition of h^^. 
First, according to the tensorial behaviour on )C", h^^ can be divided 
into scalars hab, vectors hai and a tensor hij. The vector and tensor 
components can be further decomposed as[6l l3] 

= A/ia + /iiV; 5'/iiV=o, (7) 

htj = /iL7ij + hrtj-, h^Tj = 0, (8) 

hnj = (^b,bj - ^1^J^ h^T + '2b^,h^^]^ + h%; 

&h^^j ^0, b'h^^l^O, (9) 

where Di is the covariant derivative with respect to the metric 7^-, on IC" 
and A = D ■ D. By this decomposition, we obtain the following three 
groups of variables: 

• Scalar-type variables: hab, ha, hn, hlj^\ 

• Vector-type variables: /i^^', h^Ti y 

(2) 

• Tensor-type variable: /i^/^ . 

The Einstein equations written in terms of these variables are divided into 
three subsets each of which contains only variables belonging to one of the 
above three sets of variables. Further, the tensorial indices for vector-type 
and tensor-type variables can be eliminated if we expand these variables 
in terms of harmonic tensors on IC", and gauge- invariant variables can 
be easily constructed from the harmonic expansion coefficients, as we will 
show below. In this way, the Einstein equations for perturbations of each 
type can be reduced to a set of gauge- invariant equations with a small 
and n- independent number of entries^. In particular, for tensor pertur- 
bations, this procedure gives a single wave equation for a single variable. 



3 Stability Analysis 

3.1 4D Schwarzschild black hole 

In a four-dimensional background spacetime with a Schwarzschild black 
hole, it was shown by Regge and Wheeler[7j and Zerilli 8 that by the 
harmonic expansion on /C^ = and the Fourier transformation with 
respect to time, hab ~ fab{r)e~^'^*Yl^{6,<f)), - ■ ■ , the perturbation equa- 
tions for any type of perturbations can be eventually reduced to a single 
Schrodinger-type second-order ODE 

—^+V{r)^ = Lo^^, (10) 

where r« = dr / f{r) and the master variable <1> is a combination of 
the harmonic expansion coefficients hab, ■ ■ ■ ■ The effective potential V — 
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Vv for vector perturbations (axial modes) is given by the Regge- Wheeler 
potential 

+ (11) 



and the effective potential V = Vs for scalar perturbations (polar modes) 
is given by the Zerilli potential 

^^-^^^('^ ("^ + 2) + ^— + + (12) 

where m = (Z - 1)(Z + 2)(Z = 2, 3, • • • ) and H = m + &M/r. A tensor-type 
perturbation does not exist for n = 2. 

Since the effective potentials Vv and Vs are positive outside the hori- 
zon (r > 2M), the master equations for <1> have no regular and normal- 
izable solution for u)'^ < 0. Hence, the 4D Schwarzschild black hole is 
perturbatively stable[5]. 

Similar master equations for 4D black holes with A 7^ were derived 
by Cardoso and Lemos[TCTl With the help of these equations, the 

stability of these black holes were proved by us in Ref.|5||l|. 



3.2 Vector perturbation in higher dimensions 

Now, we show that even in higher dimensions, we can find a gauge- 
invariant master variable in terms of which the perturbation equations of 
each type can be reduced to a single second-order ODE of the Schrodinger 
type. We first treat vector perturbations. 

Vector perturbations can be expanded in terms of vector harmonics 
on /C" satisfying 

{D ■ D + kl)Y, ^ 0; &V,^0, (13) 
and harmonic tensors derived from them, 

v., =-^(aV,+£),V,), (14) 

as 

hab = 0, ha, = r/a(?/)V,, h,, = 2r'^ Ht {y)Y ij . (15) 

Here, in order to guarantee the completeness of the vector harmonics, we 
have assumed that —D ■ D is extended to a self-adjoint operator in the L^- 
space of divergence-free vector fields on A3" by the Friedrichs extension. 
Note that fc^ > in general, and for if = 1, fc^ > n - ifHlO). 

For generic modes with ky ^ (n — 1)K, we can show that the combi- 
nations 

Fa := fa + ^DaHr (16) 
kv 

give a basis for gauge-invariant variables, where Da is the covariant deriva- 
tive with respect to the metric gab on the two-space ■ In contrast, for 
the exceptional modes with ky = {n — 1)K {K = 0, 1), Yij vanishes 
identically and Ht is not defined. Hence, we obtain a smaller number 
of independent gauge-invariant variables. Since we can show that these 
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Figure 1: Examples of Vv for A = and A < 0. 



exceptional modes correspond to adding rotation with some angular mo- 
mentum and are not dynamical[5], we do not discuss them in this article. 

The vacuum Einstein equations for a vector perturbation are expressed 
in terms of Fa as 

Da {r"+\'^D,{F,/r)) - mvr"-^e,5F' = 0, (17) 
kvDa{r"-^F^)=Q, (18) 

where mv = ky — {n ~ 1)K. It is easy to see that these equations are 
equivalent to 

= e"''D6 (^r"/"*) , (19) 

-'^ + Vv{t)^ = uj^-^, (20) 
art 



where 



W = ^ 



n(n + 2)K n(n - 2) , , Sn^M 
my + — 5^ -\r^ 



2r"- 



(21) 



Although the effective potential Vy is not positive definite for large 
n, as is shown in Fig^ we can show that there exists no unstable mode, 
i.e., no L^-normalizable eigenfunction "1> with o;^ < by the following 
argument, which we call the S- deformation. 

Let I be the range of r, corresponding to the regular region outside 
of the horizon, — oo < r« < r^ao- Here, r^^s is -l-co for A > 0, but is finite 
for A < 0. Then, in the space Co°(7) of smooth functions with compact 
support, the operator 

^:=-|i + nO (22) 
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is symmetric. Wc assume that A is extended to a self-adjoint operator 
in L^{I) by the Friedrichs extension. Then, the lower bound for the 
spectrum of A coincides with the lower bound of 

= $), $ e C^{I). (23) 

Here, for any regular function S{r) on I and $ € Co°{I), a partial inte- 
gration yields 



(3>, A^) = J dr, 
where 



(24) 



^=rf^+5«, (25) 

V = V + f^-S''. (26) 
ar 

Hence, if we can show that V is non-negative for an appropriate S, we 
can conclude that liP > 0. 

If we apply this method to a vector perturbation using S = n//(2r), 
we obtain 

Vv = mvL (27) 

Since mv = ky — {n — \)K is always non-negative, this implies that 

the black hole is stable against the vector perturbation in any spacetime 
dimension D = n + 2 > 4, irrespective of the values of K and A. 



3.3 Scalar perturbation in higher dimensions 

Next, wc consider scalar perturbations. Scalar perturbations can be ex- 
panded in terms of scalar harmonics on K,^ satisfying 

[a + fe^ j S = 0, (28) 

and harmonic vectors and tensors derived from them, 

Si:=~bS, (29) 

Sii - T2^'^^^ + -^iJ^' = 0' (30) 

K Ti 

as 

hab = fabS, hat=rfaS>i, hij = 2r'^ (H L'JijS + HtS^j) . (31) 

Some comments are in order. First, since we are assuming that /C" is 
complete, —A = —D ■ D has a unique self-adjoint extension in L^(/C"), 
which coincides with the Friedrichs extension and is non-negative. Second, 
we do not consider the zero modes (fc^ = 0), by assuming that such a mode 

corresponds to a variation of the parameters of the background solution. 
This assumption is satisfied when /C" is compact and closed. Third, for 
K = 1, the second smallest eigenvalue is = n. For this eigenvalue, Sij 
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vanishes. Since we can show that such modes are pure gaugef^, we do 
not consider them in this article. 

From gauge transformation laws of the perturbation variables, we find 
that the following combinations can be used as a basis for gauge-invariant 
variables pQ: 

F~Hl + -Ht + —X\ Fab ~ fab + DaXb + DbXa, (32) 

n r 

where 

Xa ■- ^ (/a + ^DaHr) ■ (33) 
In terms of the gauge-invariant variable <E> constructed from these, 

*.I!!^A^,I£). (34) 



H \ iior ^ 

n(n+l)M ,2 
H = ra+— — , m = k-nK, (35) 

yn — l ' ' ^ ' 

the Einstein equations are reduced to the single equation |5] 



' drl 



_ + Vs(r)$ = a;'<E>. (36) 



The effective potential Vs is given by 

Vs(r) = (37) 

with 

U{r) = - [n^{n + 2)(n + ifx'^ - Vln(n + l)(n - 2)mx 
+4(n - 2)(n - 4)m^] Ar^ + n(n + ifji^ 
+n{n + 1) [4(2n^ - 3n + 4)m + n(n - 2)(n - 4)(n + 1)7^] 
-12n [{n - 4)m + n{n + l)(n - 2)K] mx 

+16m^ + 4Kn{n + 2)m^ , (38) 

where x = 2M/r""\ 

We can easily see that for n — 2 and K = 1, Vs is positive definite. 
Hence, 4D spherically symmetric neutral black holes are perturbatively 
stable, irrespective of the value of the cosmological constant. In contrast, 
Vs is not positive definite for larger n in general, as illustrated in Fig|21 
Nevertheless, for K = 1 and A = 0, neutral black holes are perturbatively 
stable in any spacetime dimensions, because the S'-deformation of Vs by 

S^l^; h = r'''^+''^H (39) 
h dr 

yields 

Vs = ^ijj [ln{n + l)x + 2{l - 1) {n' n{?,l - 2) + {I ~ if}] > 0. 



(40) 
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Figure 2: Examples of Vs for A = and A < 0. 



3.4 Tensor perturbation 

In the case in which JC" is maximally symmetric, the only non-trivial ten- 
sor on /C" is the metric "fij. Hence, only the Laplacian appears as the 
derivative operator in the perturbation equation, and by the harmonic 
expansion with respect to D ■ D, the Einstein equations reduce to equa- 
tions on Af^ for any type of perturbation. In contrast, when IC" is a 
generic Einstein space, the Weyl curvature couples tensor perturbations, 
and such a reduction cannot be achieved by the harmonic expansion. How- 
ever, since the covariant derivatives and the Weyl curvature appear in the 
Einstein equations for tensor perturbations as the combination called the 
Lichnerowicz operator|12| 

ALh,j ■- -D ■ DKj - 2R^kjih''' + 2(n - l)Kh,j, (41) 

we can utilise the eigentensors of A_l, 

AlT.j = Tl = 0, D'T^j = 0. (42) 

Tensor perturbations can be expanded in terms of these eigentensors as 

hab = 0, ha^ = 0, h,j ^ 2r^ Ht {y)T ,j . (43) 

Ht is gauge-invariant by itself, and the Einstein equations are reduced 

to 

OHt + -Dr ■ DHt — Ht = 0. (44) 

In terms of 

$ = r-^/^ffr, (45) 
this equation is rewritten in the canonical form 

~~-rY + VT{r)^ = ui^^, (46) 
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where 



Vt = 



/ 



nrf , n(n-2)/ 



(47) 



Note that this equation holds for any /(r), if the background metric takes 
the form @. 

For /(r) given by 1471 can be written as 



Vt = 4 



(48) 



When K." is maximally symmetric, A_l is related to the eigenvalue fcy of 
-D ■ D hy Xl = + 2nK, and Xl - (3n - 2)K = {I - + n) > 
for K = 1. Hence, Vt > and the black hole is stable against the tensor 
perturbation. 

In contrast, when /C" is a generic Einstein space, no general lower 
bound for Xl is know, and generalised static black holes can be unstable 
against a tensor perturbation 



4 Extension to Charged Black Holes 

The extension of the formulations by Regge- Wheeler and Zerilli to Reissner- 
Nordstrom black holes was done by Moncrief, Zerilli and Chandrasekhar[T^ 
1141 115| . and used to prove their stability. In the charged black hole case, 
we have to take into account perturbations of electromagnetic fields, which 
couple metric perturbations non-trivially. However, by taking appropriate 
linear combinations of variables describing these perturbations, the per- 
turbation equations for the Einstein-Maxwell system can be reduced to 
decoupled master equations, each of which is a single second-order ODE 
of the Schrodinger equation type. Now, we show that such a reduction 
is also possible in the case of higher-dimensional static black holes with 
charge, and discuss their stabilityj^. 

We assume that an background spacetime is again locally written as 
« A/"^ X /C", and the metric is given by Q with ||HJ. We further 
assume that the electromagnetic tensor Tfi.v for the background EM field 
takes the form 

r = ^Eoeabdy" A dy\ (49) 

Then, the Maxwell equations and the Einstein equations determine f{r) 
and Eo{y) as 

fir)^K-Xy^-^ + -g^, (50) 
r" n(n — 1) 

4.1 Tensor perturbation 

Since the EM field does not couple tensor perturbations of the metric, the 
master equation for the charged case is again given by 14611 with 14711 . If 
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we apply the 5'-deformation with S — — n//(2r) to Vr, we obtain 

Vt = ^ [Al - 2{n - 1)K] . (52) 

Hence, a generalised static black hole with charge is still stable against 
tensor perturbations, if \l > 2{n — 1)K, or kj- > ~2K when /C" is 
maximally symmetric. This condition is satisfied by black holes with flat 
or spherical horizons. However, the case -K" = — 1 (A < 0) may be unstable 
even if /C" is maximally symmetric. 



4.2 Vector perturbation 

From the Maxwell equation Vf^tJ.TvA] = 0, a vector perturbation of the 
EM field is expanded in terms of the vector harmonics as 



STab = 0, 5J^a^ = DaAV^, = A [D,Y, - D,Y,j . (53) 

Hence, a vector perturbation of the EM field is described by the single 
variable A{y), which is gauge- invariant by itself. For generic modes my = 
kY — {n—l)K 7^ 0, in terms of A and Fa, the Einstein equations are written 
as 

r-D.{^LD^n)-^^^'-^A, (55) 

and the Maxwell equation 5{WvJ-'^'^) — are written as 

^D.ir^-^D^A) - 1 (fc^ + (n - 1)K + ^^^(^) A ^ S^n. 

(56) 

As shown in by taking linear combinations $± — a±r "^^Q + 
h±r''/^-^A with some constants a± and 6±, these perturbation equations 
for the Einstein-Maxwell system are transformed to the two decoupled 
equations 

+ VV±(r)$± = w^$±. (57) 



drl 

The effective potentials W± are given by 



Vv± = 4 



,2 (n^-2n + A)K n(n - 2) ^ 2 , w(5n_-2)Q^ «±_ 



(58) 

where 

^± := -!^^±1m ± A; A — [(71^ - ifM^ + 2n{n - l)mvQ^] ^'"^ . 

^ (59) 
In the limit (3 = 0, "I>+ and <&_ coincide with A and il, respectively. 
Hence, "1>+ and represent the EM mode and the gravitational mode 
of the perturbation, respectively. 
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By the 5- deformation with 5* = ^ , Vv± are transformed to 



VV± = 4 



2r • 

^2 



(n^ - 1)A/ ± A 
mv + - ' 



p-n, — 1 



(60) 



Hence, VV+ is always positive. However, since A > (n^ — 1)M for mv > 0, 
W- can become negative. 

By examining the behaviour of VV-, we find that it is positive outside 
the horizon for K = Q and K = \, while it becomes negative near the 
horizon for if = —1 if A is sufficiently close to Ac-, where Ac- is the lower 
bound on A for a fixed value of Q when the spacetime contains a regular 
black hole0. Hence, we can conclude that a charged black hole with a flat 
or spherical horizon is stable against vector perturbations. However, we 
cannot prove the stability of a charged black hole with hyperbolic horizon 
by this simple method. 

4.3 Scalar perturbation 

A scalar perturbation of the EM field is expressed as 

+ Dc(SoX'=)e„5§ = ^ {~eA + y^Ft - 2nF)) §, (61) 

&Ta^ - kEoeatX''S, = eaiD" AS^ , (62) 
STij = 0. (63) 
The expansion coefficient A is gauge invariant and obeys the equation 

r-'D. (^) - = ^^(F: - 2nF). (64) 
The Einstein equations for the scalar perturbation are equivalent to 

fD^D''^ - V5(r)$ = If^^EoA, (65) 
where Psi{r) is a function of r, and 

^^^r^ + 2i.), (66) 
H = m-\ — -; m = k ~nK. (67) 

^n — 1 ^2n — 2 ' ^ ' 

In terms of the variables $± written as linear combinations of $ and 
A, these perturbation equations can be rewritten as the two decoupled 
equations 

+ Vs±{r)^± =uj^^±. (68) 



Here, the effective potentials Vs± are given by 

1± 



- 6^5%- (69) 
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where in terms of the positive constant 6 defined by 



= {n + lfM'^6{l + mS), 



(70) 



H± are expressed as 



n{n + l)SM 



, H- =m + 



n{n + + m6)M 



(71) 



and U± are expressed in terms of a; = as 

U+ = [-4n^(n + 2)(n + ifd^x^ - 48n^(n + l)(n - 2)(5a; 
-16(n-2)(n-4)] Ar^ 
-<5^n^(3n - 2)(n + 1)''(1 + m5)x* 
+4S'^n^{n + if {{n + l)(3n - 2)m<5 + 4n^ + n - 2} 
+46{n + 1) {(n - 2)(n - 4)(ri + l)(m + n^ii:)<5 

-7n^ + 7n^ - 14n + 8} 
+ {I6(n + 1) (-4m + 3n^(n - 2)/r) 5 - 16(3n - 2)(n - 2)} x 
+64m + 16n(n + 2)K. (72) 

?7_ = [-4n^(n + 2)(n + 1)^(1 + m5)^a;^ 

+48n^(n + l)(n - 2)m(l + m5)x - 16(n - 2)(n - 4)m^] Ar^ 
-n'^(3n - 2)(n + 1)*(5(1 + m<5)^x'' 
-4n\n + + mSf {{n + l)(3n - 2)m5 -~ n^} x^ 

+4{n + 1)(1 + m5) {m(n - 2)(n - 4)(n + l)(m + n^K)5 
+4n(2n^ - 3n + 4)m + n^(n - 2)(n - 4)(n + 1)A'} a;^ 
-16m {(n + l)m (-4m + 3n^(n - 2)K) 5 
+3n{n - 4)m + 3n^(n + l)(n - 2)K} x 

+64m^ + 16n(n + 2)m^i^. (73) 
By the S'-deformation with 



Since this is positive definite, the EM mode is always stable. 
By the similar S-deformation 



(74) 



the effective potential Vs+ for is transformed to 



Vs+ = ,^^\{n-2){n + \)5x + 2]. 



(75) 



the effective potential Vs- for is transformed to 




(76) 



^s- = TAir- [2rn - (n + l)(n - 2)(1 + m6)x] . 



(77) 
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Figure 3: Examples of Vs- for A = and A < 0. 



Table 1: Stability of static black holes. 
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D > 6 ? 
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> 


OK 


OK 


OK 


OK 


D < 6 OK 
D > 7 ? 


-D = 4, 5 OK 
D > G 7 
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< 


OK 


OK 


OK 


OK 


i) = 4 OK 

D > 5 ? 


L> = 4 OK 
r> > 5 ? 


K = 


X 


< 


OK 


OK 


OK 


OK 


D = 4 OK 
D > 5 ? 


D = 4 OK 
r» > 5 ? 


K = -1 


X 


< 


OK 




OK 




n = 4 OK 
D > 5 ? 


O = 4 OK 
D > 5 ? 



Hence, the gravitational mode is also stable for n — 2. 

However, for n > 2, Vs- becomes negative near the horizon in general 
(See FigOJ. The S'-deformat ion used in the neutral case does not work 
either, because it cannot remove a negative ditch of Vs- produced by the 
charge near the horizon. We have not been able to prove the stability for 
n > 2 by this method, except for some special cases. 

5 Summary 

In this article, we have discussed the stability of generalised static black 
holes with and without charge in higher dimensions. The results of our 
analysis for maximally symmetric black holes are summarised in Table^ , 
In particular for the case of A = 0, we have established the stability of 
a spherically symmetric black hole for arbitrary dimensions, and that of 
of a spherically symmetric black hole with charge for D = 4 and D = 5. 
We have also proved the stability of 4D and 5D dS-Schwarzschild and 
4D AdS-Schwarzschild black holes with charge. Although we have not 

^In this table, the results for vector and scalar perturbations are also valid for generalised 
static black holes, except for the case K = I, \ > and Q = 0, for which the stability has 
been proved for _D < 5 in spherically symmetric case. 
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succeeded in proving the stability of black holes for the other cases, there 
is no strong reason to suspect the instability of such black holes because 
the negative ditch of the effective potential is not so deep. 

Finally, we would like to point out that the master equations used in 
this article for the stability analysis can be also used in other problems 
related to black holes, such as perturbative analyses of the uniqueness 
of asymptotically de Sitter or anti-de Sitter static black holes in four 
and higher dimensions, the precise determination of the frequencies of 
quasinormal modes and the greybody factor of the Hawking radiation, 
and the estimation of the gravitational wave emission rate from mini black 
holes, which might be produced in colliders and cosmic ray events^. 
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